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Abstract. We initiated a comprehensive state of the art binary population synthesis study of white dwarf main- 
sequence star (WDMS) binaries to serve as a foundation for subsequent studies on pre-cataclysmic variables, double 
white dwarfs, and white dwarf + B-star binaries. We considered seven distinct formation channels subdivided into 
three main groups according to the evolutionary process that gives rise to the formation of the white dwarf or 
its helium-star progenitor: dynamically stable Roche- lobe overflow (Algol- type evolution), dynamically unstable 
Roche- lobe overflow (common-envelope evolution), or stellar winds (single star evolution). For each formation 
channel, we examine the sensitivity of the population to changes in the amount of mass lost from the system 
during dynamically stable Roche-lobe overflow, the common-envelope ejection efficiency, and the initial mass 

■ ratio or initial secondary mass distribution. In the case of a flat initial mass ratio distribution, the local space 
density of WDMS binaries is of the order of ~ 1CP 3 pc -3 . This number decreases to ~ 10~ 4 pc~ 3 when the initial 
mass ratio distribution is approximately proportional to the inverse of the initial mass ratio. More than 75% of 
the WDMS binary population originates from wide systems in which both components essentially evolve as if 
they were single stars. The remaining part of the population is dominated by systems in which the white dwarf 
is formed in a common-envelope phase when the primary ascends the first giant branch or the asymptotic giant 
branch. When dynamically stable mass transfer proceeds highly conservative and the common-envelope ejection 
process is very efficient, the birthrate of WDMS binaries forming through a common-envelope phase is about 10 
times larger than the birthrate of WDMS binaries forming through a stable Roche-lobe overflow phase. The ratio 
of the number of helium white dwarf systems to the number of carbon/oxygen or oxygen/neon/magnesium white 
dwarf systems derived from large samples of observed WDMS binaries by, e.g., future planet-search missions such 
as SuperWASP, COROT, and Kepler may furthermore constrain the common-envelope ejection efficiency. 
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1. Introduction 0.3 M are generally harder to detect due to the large lu- 
minosity difference between the component stars (Iben et 

In recent years, substantial numbers of detached white al 1997; Margh 2000 )_ Despite this obstacle, several claims 

dwarf main-sequence star (WDMS) binaries have been de- have recently been made on the possible detection of white 

tected in large-scale surveys searching for pre-cataclysmic dwarf binaries with bright companions of spectral types as 

variables (Hillwig et al. 2000, Raymond et al. 2003), dou- early as B(e) (e g _ Ho lberg et al. 1998, Vennes et al. 1998, 

ble degenerates (Saffer et al. 1998), low-mass white dwarfs Burleigh & Ba rstow 2000, Burleigh et al. 2001, Shobbrook 

(Marsh et al. 1995, Maxted et al. 2000), planetary nebulae e ^ a j 2003' and references therein) 
(Bond & Livio 1990, Livio 1992), or even the dark matter 

content of the Galaxy's heavy halo (Silvestri et al. 2002). WDMS binaries with early B(e)-type companions 

Upcoming surveys looking for transiting extrasolar plan- more massive than ~ 10 M Q are thought to be the end 

ets are furthermore expected to contribute further to this product of binaries in which a giant-type star transfers 

rapidly growing database (Farmer & Agol 2003). So far, mass to its initially less massive companion. The mass- 

the majority of WDMS binaries found consist of a white transfer phase is responsible for exposing the giant's core 

dwarf with a low-mass companion. White dwarfs in bi- as a white dwarf and for increasing the mass of the sec- 

naries with main-sequence secondaries more massive than ondary so that it evolves towards earlier spectral types. In 

addition, transport of angular momentum may spin the 

Send offprint requests to: B. Willems secondary up to form a rapidly rotating Be star (Waters 

* Present address: Northwestern University, Department of et al. 1988, Pols et al. 1991). The predicted number of 

Physics and Astronomy, 2145 Sheridan Road, Evanston, IL WDMS binaries forming in this way is, however, rather 

60208, USA small (e.g. de Kool & Ritter 1993, Figs. 2 and 4) which 
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may pose a problem if the recent claims on the detection 
of white dwarf + B-star binaries are confirmed. In addi- 
tion, as mass transfer is expected to increase the orbital 
period it is hard to explain the formation of short-period 
white dwarf + early B-star candidates such as HD 161103 
(-forb = 4.7 days, Shobbrook et al. 2003) through this evo- 
lutionary channel. 

The wealth of existing and expected new data com- 
bined with new up-to-date stellar evolution tools (see 
Hurley et al. 2000) render the WDMS binaries a poten- 
tially fruitful subject of a comprehensive binary popula- 
tion synthesis study. A systematic exploration of the for- 
mation channels in particular may shed new insights on 
the origin of short-period white dwarf + B-star binaries. 
Other interesting subclasses of detached WDMS binaries 
include pre-cataclysmic variables and progenitors of dou- 
ble degenerates. 

To the best of our knowledge, de Kool & Ritter (1993) 
were the first to devote an entire study to the population 
of WDMS binaries as a class on its own. The authors used 
a Monte Carlo type simulation in which the initial binary 
parameters were randomly generated from observed or- 
bital period, primary mass, and secondary mass distribu- 
tions; and the evolution of the binaries was approximated 
by analytical fits to detailed stellar evolution tracks as 
described by Politano (1988) and dc Kool (1992). The 
method resulted in distribution functions for the forma- 
tion rates and space densities of detached WDMS binaries 
(or their direct progenitors) as a function of their expected 
orbital and stellar parameters. 

In this paper, our aim is to update the population 
synthesis study by de Kool & Ritter (1993) by using the 
more recent and much more detailed analytical fits to stel- 
lar evolution derived by Hurley et al. (2000). We study 
in more detail the differences between the evolutionary 
channels leading to the formation of WDMS binaries, and 
derive the formation rates and absolute numbers of sys- 
tems for each of the formation channels separately. In ad- 
dition, we consider a much more elaborate set of popula- 
tion synthesis models characterised by different assump- 
tions about the fate of the mass transferred during dy- 
namically stable Roche-lobe overflow, different common- 
envelope ejection efficiencies, and different initial mass ra- 
tio or secondary mass distributions. 

Our method furthermore differs from that by de Kool 
& Ritter (1993) in the implementation of the observed dis- 
tributions of the initial masses and orbital periods: instead 
of randomly generating them from the observed distribu- 
tion functions at the onset of the binary evolution calcu- 
lations, we start our simulations from a logarithmically 
spaced 3-dimensional grid of initial masses and orbital pe- 
riods in which each set of binary parameters has an equal 
probability. The likelihood of the different initial binary 
parameters is then taken into account by weighting the 
contribution of each WDMS binary found in the simula- 
tions according to the adopted primary mass distribution, 
initial mass ratio or secondary mass distribution, and ini- 
tial orbital period distribution. In doing so, we obtain a 



better and more uniform sampling of the entire parameter 
space, which particularly benefits the resolution of low- 
density tails in the distribution functions describing the 
masses and orbital periods of WDMS binary components. 
This procedure also allows us to vary the distribution func- 
tions for the initial masses and orbital periods without 
having to repeat the binary evolution calculations. 

The results of this study are intended to serve as a 
foundation for more focused population synthesis studies 
of white dwarf + B-star binaries (Willems et al. , in prepa- 
ration), double white dwarfs (Willems et al., in prepara- 
tion), and pre-cataclysmic variables. Hence, we here con- 
centrate on the effect and the uncertainties associated 
with the different model parameters in general. In sub- 
sequent papers we hope to constrain these uncertainties 
by comparing our results with different types of WDMS 
binaries and WDMS binary descendants. The results pre- 
sented here therefore do not necessarily represent the best 
possible fit to observationally derived binary parameters. 

The plan of the paper is as follows. In Section 2, 
we present the basic ingredients and the assumptions 
adopted in the binary population synthesis calculations. 
In Section 3, we describe the different evolutionary chan- 
nels leading to the formation of WDMS binaries as well as 
the ranges of stellar masses and orbital periods available 
to each formation channel. Section 4 deals with the ef- 
fects of varying degrees of non-conservative mass transfer 
and different common-envelope ejection efficiencies on the 
population of WDMS binaries. In Section 5, we estimate 
the order of magnitude of the formation rates and the ex- 
pected number of WDMS binaries currently populating 
the Galaxy. Section 6 briefly addresses the expected lumi- 
nosity differences between the white dwarf and its com- 
panion, as well as their expected orbital radial-velocity 
variations. A brief summary of our results and some con- 
cluding remarks are presented in Section 7. 

2. Basic concepts and assumptions 

2.1. The population synthesis code 

We use the BiSEPS binary population synthesis code in- 
troduced by Willems & Kolb (2002) to study the forma- 
tion of detached white dwarf main-sequence star binaries. 
The code uses the single star evolution formulae derived 
by Hurley et al. (2000) and follows the main steps of the 
binary evolution scheme outlined by Hurley et al. (2002). 
All binary orbits are treated as circular and stellar rota- 
tion rates are kept synchronised with the orbital motion at 
all times. For the purpose of this investigation, we further- 
more limit ourselves to Population I stellar compositions. 

When a binary component overflows its Roche-lobe, 
the stability of the ensuing mass-transfer phase is deter- 
mined by means of the radius-mass exponents introduced 
by Webbink (1985). Mass transfer taking place on the dy- 
namical time scale of the donor star is assumed to lead to 
a common-envelope (CE) phase during which the orbital 
separation is reduced and the envelope is expelled from 
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the system. The phase is modelled in the usual way by 
equating the binding energy of the donor star's envelope 
to the change in the orbital energy of the binary compo- 
nents. The orbital separation af at the end of the common- 
envelope phase is then related to the orbital separation a; 
at the beginning of the phase as 



Of 

a; 



M c /M 1 



1 + 2 M e / (M2 «ce Atl) ' 



(1) 



Here G is the gravitational constant, M c and M e are the 
core and envelope mass of the Roche- lobe overflowing star, 
Mi = M c + M c is the total mass of the donor star, M2 
is the mass of the companion, and tl = Rl/cli is the ra- 
dius of the donor star's Roche lobe in units of the orbital 
separation at the start of the common-envelope phase. For 
the binding-energy parameter A and the common-envelope 
ejection efficiency acE, we adopt the commonly used val- 
ues A = 0.5 and q;ce = 1-0 (e.g. de Kool 1992, Politano 
1996). 

In the case of dynamically stable mass transfer, a frac- 
tion 1 — 7rlof of the transferred mass is assumed to be 
accreted by the donor star's companion, while the remain- 
ing fraction 7rlof is lost from the system carrying away 
the specific orbital angular momentum of the companion. 
The amount of angular momentum carried away by mat- 
ter leaving the system during non-conservative Roche-lobe 
overflow is, however, still an unresolved issue, so that this 
quantity is effectively a free parameter (see, e.g., the ap- 
pendix in Kolb et al. 2001). 

For non-degenerate accretors, we set 



1 - 7RLOF = min 10 



M 



1 



THK,. 



(2) 



where is the mass-transfer time scale of the donor 
and thk, a is the thermal time scale of the accretor. With 
this prescription, the mass-transfer phase is conservative 
as long as THK,a < 10 r A ^ (see also, for example, Iben 
& Tutukov 1987, Pols et al. 1991, Hurley et al. 2002). 
Accretion onto white dwarfs is assumed to be fully non- 
conservative so that 7r,lof = 1 and M W d = 0, where 
Mwd is the mass of the white dwarf. This assumption 
does not affect the formation of detached WDMS bina- 
ries, but may influence the fate of their descendants when 
the main-sequence star becomes larger than its critical 
Roche lobe. We do not deal with neutron star and black 
hole accretors in this investigation. 

For more details on the treatment of mass-loss and 
mass-accretion in the BiSEPS code we refer to Willems & 
Kolb (2002). 

2.2. Initial masses and orbital periods 

We start our population synthesis study by evolving a 
large number of binaries initially consisting of two zero-age 
main-sequence stars with a mass between 0.1 and 30 M Q 
and an orbital period between 0.1 and 100 000 days. The 
initial primary and secondary masses, Mi and M 2 , and the 



initial orbital periods P rb are taken from a grid consisting 
of 60 logarithmically spaced stellar masses and 300 loga- 
rithmically spaced orbital periods. The maximum evolu- 
tionary age considered for each binary is 10 Gyr. For sym- 
metry reasons only binaries with Mi > Mi are evolved. 

The number of systems following a sequence of evolu- 
tionary phases similar to those of a given binary in our 
simulated sample is determined by the probability of the 
binary's initial parameters, by the star formation rate at 
the birth of the zero-age main-sequence binary, and by the 
fraction of stars in binaries. We assume the initial primary 
masses to be distributed according to the normalised ini- 
tial mass function (IMF) 1 

Mx/Mq < 0.1, 

£ (Mi) = \ 0.38415 Mf 1 0.1 < M x /M < 0.75, (3) 
0.23556 M-f 2,7 0.75 < Mi/M < 00, 

the initial mass ratios q = M2 jM\ according to 



n(q) 



[iq v < g < 1, 



q>l, 
and the initial orbital separations a according to 



(4) 



X(a) 







a/R e < 3 or a/R Q > 10 6 , 



0.078636 a- 1 3 < a/R® < 10 6 



(5) 



In Eq. v is a constant and fi a normalisation factor 
depending on v. Unless stated otherwise, we set v = and 
H = 1. The upper limit of 10 6 Rq in the distribution of 
the initial orbital separations is larger than in our previ- 
ous investigations in order to properly take into account 
the contribution of very wide systems to the population 
of WDMS binaries. For more details and references on 
the adopted distribution functions, we refer to Willems & 
Kolb (2002). 

We furthermore assume all stars to be in binaries and 
adopt a constant star-formation rate S calibrated so that 
one binary with Mi > 0.8 M is born in the Galaxy each 
year (see also Iben & Tutukov 1984, Han et al. 1995, 
Hurley et al. 2002). When combined with an effective 
Galactic volume of 5 x 10 11 pc 3 , this yields an average local 
birthrate of Galactic white dwarfs of 2 x 10~ 12 pc -3 yr _1 , 
which is consistent with observations (Weidemann 1990). 
From this calibration, it follows that 



S 



£ > {M 1 )dM 1 = 1, 



(6) 



so that S — 4.9 yr 1 . This rate may be converted into 
an approximate local star formation rate (expressed in 

We note that the adopted IMF is a simplified version of the 
IMF by Kroupa et al. (1993). The simplification is introduced 
because of the still existing uncertainties in the IMF for low- 
mass stars (see, e.g., Scalo 1998, Kroupa 2001). Its effect on 
our results is small in comparison to the overall uncertainties 
of the population synthesis models. 
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pc~ 3 yr -1 ) by dividing it by 5 x 10 11 pc 3 . As the star for- 
mation rate may have been higher in the past, the cali- 
bration of the rate to match the observationally inferred 
current birthrate of Galactic white dwarfs may yield an 
underestimate of the number of binaries with old compo- 
nent stars (Boissier & Prantzos 1999). 



3. Formation channels 

Stars in close binaries can evolve into white dwarfs either 
through the loss of their envelope by the action of a stel- 
lar wind or by mass transfer resulting from dynamically 
stable or unstable Roche-lobe overflow. For brevity, we re- 
fer to binaries that do not undergo mass-transfer episodes 
as non- interacting binaries, even if some mass exchange 
and orbital evolution takes place due to the action of a 
stellar wind. Binaries in which a white dwarf or its direct 
progenitor is formed as the end product of mass trans- 
fer are referred to as interacting binaries. In what follows, 
we divide the latter group according to the stability of 
the mass-transfer phase and according to the remnant it 
leaves behind. 

Since detached WDMS binaries constitute an impor- 
tant intermediary phase in the formation of many more 
exotic binary systems such as cataclysmic variables and 
double degenerates, some of the formation channels de- 
scribed below have been partially discussed before, albeit 
maybe in lesser detail and in different contexts. The evolu- 
tion of binaries with low- to intermediate-mass component 
stars has been studied in detail by, e.g., Iben & Tutukov 
(1985, 1986a, 1987), van der Linden (1987), de Loore & 
Vanbeveren (1995), Langer et al. (2000), Han et al. (2000), 
Nelson & Eggleton (2001), Chen & Han (2002, 2003); and 
references therein. 

3.1. Dynamically stable mass transfer 

3.1.1. Case B RLOF with a white dwarf remnant 

The first formation channel applies to initial binaries con- 
sisting of two low-mass main-sequence stars with orbital 
periods that are too short to allow the primary to evolve 
on the giant branch without overflowing its Roche lobe. 
Most of the systems initiate mass transfer on the thermal 
time scale of the donor star when it approaches the end of 
the main sequence or when it crosses the Hertzsprung gap. 
Once the mass ratio is inverted, mass transfer generally 
slows down and continues as the donor star ascends the 
first giant branch. In most cases, the case B mass-transfer 
phase is more important for determining the subsequent 
evolution than the initial case A phase. For brevity, we 
therefore simply refer to the case B phase as the phase 
characterising the formation channel. 

Since the thermal time scale of the accreting com- 
panion is of the same order of magnitude as that of the 
donor star, the mass-transfer phase is highly conservative 
[see Eq. (j2J]. The secondary's mass and the orbital pe- 
riod may therefore increase substantially with respect to 
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Fig. 1. Schematic representation of the main evolutionary 
phases leading to the formation of a WDMS binary via 
evolutionary channel 1 (see text for details). MS stands 
for main sequence, HG for Hertzsprung gap, GB for giant 
branch, and WD for white dwarf. 



their values at the onset of the Roche-lobe overflow phase. 
The mass-transfer phase terminates when the giant's en- 
tire envelope is transferred to the companion, exposing its 
helium core as a low-mass white dwarf. The main evolu- 
tionary phases characterising this formation channel are 
summarised schematically in Fig. ^ For future reference, 
we label this channel as channel 1. 

The two-dimensional probability distribution func- 
tions (PDFs) describing the regions of the (Mi,P or b)- 
and the (Ma, P or b)-planes occupied by binaries evolving 
through formation channel 1 are displayed in the upper- 
panels of Figs. [3 and The left-hand panels of the fig- 
ures show the distributions of the WDMS binary progeni- 
tors at the beginning of their evolution as zero-age main- 
sequence binaries, while the right-hand panels show the 
distributions of the WDMS binaries at the time of their 
formation. For the construction of the plots, only bina- 
ries contributing to the present-day Galactic population 
are taken into account 2 . The distributions are normalised 
so that the bin containing the largest contribution to the 
population has a PDF value equal to one. In order to 
show the maximum amount of detail this normalisation 



2 Note that our purpose here is to examine the formation 
space of WDMS binaries. The plots are therefore not fully 
representative for the present-day population of short-period 
systems (P rb 1 day) which may have undergone significant 
orbital shrinkage due to magnetic braking and/or gravitational 
radiation after their formation. However, even though the or- 
bital evolution is not shown in the figures, the finite life time 
of the systems resulting from the orbital evolution and/or the 
nuclear evolution of the secondary is taken into account in the 
determination of the WDMS binaries currently populating the 
Galaxy. 
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Fig. 2. Distribution of present-day WDMS binaries in the (Mi, P or b)-plane at the beginning of the progenitors' evo- 
lution (left-hand panels) and at the beginning of the WDMS binary phase (right-hand panels). The distributions are 
normalised so that for each formation channel the bin containing the largest contribution to the population has a PDF 
value equal to one. 



6 



B. Willems and U. Kolb: Detached white dwarf main-sequence star binaries 



is performed separately for each of the formation chan- 
nels considered. The relative importance of the different 
formation channels will be illustrated in Section 13.41 and 
addressed in more detail in Sectional 

The initial orbital periods of the binaries evolving 
through formation channel 1 typically range from 0.5 to 
3 days. Binaries with shorter orbital periods tend to con- 
verge before the primary is able to evolve into a white 
dwarf due to the angular momentum losses caused by mag- 
netic braking and/or gravitational radiation (Pylyser & 
Savonije 1988, 1989). Systems with longer orbital periods 
on the other hand have somewhat too evolved donor stars 
to initiate a dynamically stable Roche-lobe overflow phase. 
Instead, they undergo a common envelope phase which, in 
view of the short initial periods, generally results in the 
merger of the donor star's core with its main-sequence 
companion. 

The initial primary and secondary masses available to 
formation channel 1 are limited to the intervals given by 
1 M Q <; Mi & 5 M e and M 2 <; 3 M . The primary mass 
interval arises from the dual requirement that the primary 
must be massive enough to evolve away from the zero-age 
main sequence within the imposed age limit of 10 Gyr, but 
not so massive that its core mass at the end of the mass- 
transfer phase is high enough to ignite helium in its central 
layers. In our model, stars more massive than ~ 2.5 Mq 
are still able to evolve into helium white dwarfs due to the 
reduction in mass caused by the thermal time scale mass- 
transfer phase on the main sequence. The upper limit of 
3 Mq on the secondary mass is imposed by the secondary's 
main-sequence life time which needs to be long enough to 
allow the primary to evolve into a white dwarf before the 
secondary leaves the main sequence. Since the life time 
of a star on the main sequence decreases with increasing 
values of its mass, conservative mass transfer here yields a 
competitive race between the formation of the white dwarf 
and the accelerating evolution of the secondary. A similar 
behaviour can be deduced from the evolutionary scenarios 
described by Iben & Tutukov (1985). 

The formation channel eventually gives rise to WDMS 
binaries consisting of a 0.1 — 0.4M a He white dwarf and 
a main-sequence star with a mass up to 6 Mq . Due to 
the stable mass-transfer phase on the giant branch the or- 
bital periods may be as long as 100 days. The majority 
of the newly formed WDMS binaries have an orbital pe- 
riod which is correlated with the mass of the white dwarf. 
The correlation arises during the mass-transfer phase on 
the giant branch where the radius of the giant, which is 
approximately equal to the radius of its Roche-lobe, is de- 
termined by the mass of its core. The same relation gives 
rise to the well-known correlation between the white dwarf 
mass and the orbital period in wide binary millisecond pul- 
sars (see, e.g., Joss et al. 1987, Savonije 1987, Rappaport 
et al. 1995, Ritter 1999, Tauris & Savonije 1999). The 
small number of WDMS binaries occupying the region be- 
low the P or b — Mwd relation correspond to thermal time 
scale mass-transfer systems for which the primary already 
lost most its envelope prior to the stable mass-transfer 



phase on the giant branch. The orbital periods of the 
newly formed WDMS binaries furthermore increase with 
increasing mass of the secondary. This relation arises from 
the narrow range of initial orbital periods and secondary 
masses available to the formation channel and from the 
increase of the orbital period with the amount of mass 
transferred during the conservative Roche-lobe overflow 
phase on the giant branch. 

3.1.2. Case B RLOF with a naked helium star remnant 

Binaries evolving through the second formation channel 
start their evolution with more massive component stars 
and somewhat longer orbital periods than those in chan- 
nel 1. They again undergo a highly conservative case B 
mass-transfer phase, possibly preceded by a short case A 
phase during the final stages of the primary's evolution 
on the main sequence. Since the initial primary masses 
are higher and the initial orbital periods are longer than 
in channel 1, the core of the Roche- lobe filling giant is now 
massive enough to ignite helium in its central layers so that 
a naked helium star is formed instead of a white dwarf. 
During the subsequent evolution, the naked helium star 
at some point loses its envelope either through the action 
of a stellar wind or through a short second phase of stable 
Roche-lobe overflow. The outcome in both cases is a rela- 
tively wide binary consisting of a C/O or O/Ne/Mg white 
dwarf and an intermediate- to high-mass main-sequence 
secondary. The main evolutionary phases of this forma- 
tion channel are summarised schematically in Fig. 21 We 
will refer to the channel as channel 2. 

From the middle panels of Figs. |21 and [31 it follows 
that most of the WDMS binaries evolving through forma- 
tion channel 2 originate from binaries with initial orbital 
periods between 2 and 130 days, initial primary masses 
between 2 Mq and 12M Q , and initial secondary masses 
smaller than 11M Q . We note that there is a small num- 
ber of systems with primary masses up to 20 Mq for which 
highly conservative mass transfer results in the formation 
of secondaries with masses up to 30 Mq. However, since it 
is uncertain whether or not a white dwarf may be formed 
from primaries with such high initial masses, we do not 
include these systems in the following discussion. 

The limits on the range of initial orbital periods result 
from the evolutionary stage of the primary at the onset of 
the first Roche-lobe overflow phase. For binaries with ini- 
tial orbital periods shorter than 2 days, the core of the pri- 
mary emerges from the mass-transfer phase as a low-mass 
naked helium star with a helium-burning life time that is 
too long for it to form a white dwarf before the secondary 
leaves the main sequence. Binaries with initial orbital pe- 
riods longer than 130 days on the other hand have signifi- 
cantly evolved donor stars with deep convective envelopes. 
They are therefore subjected to a common-envelope phase 
instead of to a thermal time scale mass-transfer phase. 

The lower limit of 2 Mq on the mass of the primary 
stems from the requirement that the star must be able to 
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Fig. 3. Distribution of present-day WDMS binaries in the (M2, -P rb)-plane at the beginning of the progenitors' evo- 
lution (left-hand panels) and at the beginning of the WDMS binary phase (right-hand panels). The distributions are 
normalised in the same way as in Fig. [2] 
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Fig. 4. As Fig.Ql but for formation channel 2. nHe stands 
for naked helium star. 

develop a sufficiently massive core capable of igniting he- 
lium at the end of the stable mass-transfer phase on the 
giant branch. The upper limit of 12 Mq corresponds to 
the highest mass for which a star subjected to mass loss 
may evolve into a white dwarf rather than into a neutron 
star (e.g. van den Heuvel 1981, Law & Ritter 1983). The 
lower and upper limits on the secondary mass result from 
the requirement that mass transfer from the primary is 
dynamically stable and from our convention that the sec- 
ondary is initially less massive than the primary. 

At the time of formation, the WDMS binaries forming 
through formation channel 2 typically consist of a 0.65 — 
1.44 M Q C/O or O/Ne/Mg white dwarf and a 1 - 20 M Q 
main-sequence star. The orbital periods range from 10 to 
1000 days, with the bulk of the systems occupying a rather 
narrow band of orbital periods around P rb ~ 200 days. 

3.1.3. Case C RLOF with a white dwarf remnant 

The third evolutionary channel applies to binaries with 
initial orbits wide enough to allow the primary to evolve 
to the bottom of the asymptotic giant branch (AGB) with- 
out overflowing its Roche lobe. As it ascends the AGB, a 
strong stellar wind decreases the primary's mass below 
2/3 of the mass of its companion so that when it finally 
does fill its Roche-lobe, the ensuing mass-transfer phase is 
dynamically stable (e.g. Webbink et al. 1983). When mass 
transfer ends, any remaining surface layers are quickly re- 
moved by the wind, leading to the exposure of the AGB 
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Fig. 5. As Fig. but for formation channel 3. CHcB 
stands for core helium burning on the horizontal branch. 



star's C/O or O/Ne/Mg core as a white dwarf. The orbital 
period of the newly formed WDMS binary may be sub- 
stantially longer than the initial binary period due to the 
combined effect of the stellar wind and the stable mass- 
transfer phase on the AGB. The main evolutionary phases 
followed by a binary evolving through this formation chan- 
nel are summarised schematically in Fig. [SJ We will refer 
to the channel as formation channel 3. 

The middle panels of Figs. and |21 show that the bi- 
naries evolving through formation channel 3 have initial 
orbital periods between 600 and 10 000 days. In binaries 
with shorter orbital periods the AGB wind has insuffi- 
cient time to decrease the mass ratio below 2 /3 before the 
primary fills its Roche lobe. The resulting mass-transfer 
phase is then dynamically unstable and leads to the for- 
mation of a common envelope (see Section f3.2.3|) . Binaries 
with initial orbital periods longer than 10 000 days are too 
wide to interact. 

The initial primary and secondary masses range from 
1-M© to 8M©, similar to the mass range leading to the 
formation of white dwarfs by single stars. This conformity 
arises because the primary's Roche-lobe overflow phase 
does not take place until the very late stages of the AGB 
(see also Iben & Tutukov 1985). The fact that the mass 
ratio needs to drop below 2/3 furthermore implies that 
the initial secondary mass generally cannot be too much 
smaller than the initial primary mass. 

The formation channel eventually gives rise to WDMS 
binaries consisting of a 0.5 - 1.44 Mq C/O or O/Ne/Mg 
white dwarf and a 1 — 8 Mq main-sequence star revolving 
around each other with a period of 2000 to 20000 days. 
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Fig. 6. As Fig. ^ but for formation channel 4. 



3.2. Dynamically unstable mass transfer 

3.2.1. Case B CE phase with a white dwarf remnant 

The fourth evolutionary channel is characterised by a dy- 
namically unstable mass-transfer phase from a low-mass 
giant-branch star. As the secondary plunges into the donor 
star's rapidly expanding envelope, the orbit shrinks and 
orbital energy is transferred to the envelope until it is 
expelled from the system. At the end of the phase, the 
core of the Roche-lobe filling giant is exposed as a helium 
white dwarf which orbits the main-sequence secondary 
with a drastically reduced orbital separation. The main 
evolutionary phases comprising this formation channel are 
summarised schematically in Fig. [fj] We will refer to the 
channel as channel 4. 

It follows from the middle panels of Figs. |21 and [3] that 
the systems evolving through formation channel 4 have 
initial orbital periods between 30 and 1000 days. Binaries 
with initial orbital periods shorter than 30 days do not 
survive the common-envelope phase evoked by the pri- 
mary, while binaries with orbital periods longer than 1000 
days do not undergo Roche-lobe overflow until the pri- 
mary reaches the AGB. The limits on the initial primary 
and secondary mass ranges, 1 M <J Mi < 2 M Q and 
M 2 Si 2 M , arise for similar reasons as those in forma- 
tion channel 1 (see Section I3.1.1|) . The upper limits on 
Mi and M2 are here somewhat smaller because the longer 
initial orbital periods allow a primary of a given mass to 
reach a more evolved evolutionary state and thus to de- 
velop a more massive helium core than the same primary 
in a shorter-period binary. The mass and period ranges 
of binaries surviving the common-envelope phase in this 
evolutionary channel are in excellent agreement with the 
three-dimensional hydrodynamical simulations performed 
by Sandquist et al. (2000). 

The formation channel results in WDMS binaries con- 
sisting of a 0.3 — 0.5 Mq He white dwarf and a main- 
sequence star with a mass up to 2 M Q . The orbital period 
takes values between 0.05 and 30 days. The majority of 
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Fig. 7. As Fig. ^ but for formation channel 5. 



the systems have P or b ^ 3 days, M\vd ~ 0.4 M Q and 
M 2 £1M . 

3.2.2. Case B CE phase with a naked helium star 
remnant 

In the fifth evolutionary channel, a common-envelope 
phase occurs when an intermediate- to high-mass primary 
crosses the Hertzsprung-gap or ascends the first giant 
branch. The core of the primary emerges from the spiral- 
in phase as a low-mass naked helium star which orbits the 
main-sequence secondary with an orbital period of the or- 
der of a few days. The following few million years, the 
primary burns helium in its core until it loses its surface 
layers in a stellar wind or until it overflows its Roche lobe 
a second time when helium is exhausted in the core. In 
either case, the primary evolves into a C/O or O/Ne/Mg 
white dwarf orbiting an intermediate- to high-mass main- 
sequence secondary. The main evolutionary phases charac- 
terising this formation channel arc summarised schemati- 
cally in Fig. We will refer to this channel as formation 
channel 5. 

The upper panels of Figs.Eland|21show that the forma- 
tion channel applies to binaries with initial orbital periods 
between 10 and 1000 days, initial primary masses between 
4M© and 12 Mq, and initial secondary masses between 
1 M© and HMq. Stars in binaries with orbital periods 
shorter then 10 days are so close that Roche-lobe over- 
flow from the primary usually already occurs before it has 
developed a deep convective envelope, so that mass trans- 
fer tends to be dynamically stable (see formation chan- 
nel 2, Section f3.1.2|l . Binaries with periods longer than 
1000 days on the other hand allow the primary to evolve 
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beyond the first giant branch without filling its Roche 
lobe. The gap in the initial period range around ~ 100 
days separates systems initiating mass transfer in the 
Hcrtzsprung gap from those initiating mass transfer on the 
giant branch. Mass transfer from Hertzsprung-gap donor 
stars often starts of as a thermal time scale mass-transfer 
phase which evolves into a common-envelope phase as 
the donor star approaches the giant branch. The sur- 
vival of these systems depends strongly on the ability of 
the thermal time scale mass-transfer phase to decrease 
the primary's mass sufficiently before the onset of the 
common-envelope phase. The gap between systems with 
Hcrtzsprung-gap and giant-branch donor stars is related 
to the behaviour of the adiabatic radius-mass exponents 
tabulated by Hjellming (1989). The large values found by 
Hjcllming (1989) near the transition phase where the star 
starts to develop a deep convective envelope yield a small 
window in the parameter space where the mass-transfer 
phase is dynamically stable, so that a different evolution- 
ary scenario ensues (e.g. formation channel 2). However, in 
view of the still existing uncertainties in the detailed mod- 
elling of this transition phase, the associated values of the 
adiabatic radius-mass exponents are also quite uncertain. 
The occurrence of the gap may therefore be an artifact of 
the stability criterion separating systems undergoing dy- 
namically stable Roche-lobe overflow from those undergo- 
ing dynamically unstable Roche-lobe overflow. The origin 
of the limits on the primary mass range is similar to that of 
the limits found for formation channel 2 (Section 13.1.2(1 . 
The lower limit of 1 Mq on the secondary mass corre- 
sponds to the smallest companion mass for which the bi- 
nary is able to avoid a merger. 

At the time of formation, the WDMS binaries forming 
through formation channel 5 consist of a 0.65 — 1.44 Mq 
C/O or O/Ne/Mg white dwarf and a 2 - 13 Mq main- 
sequence star orbiting each other with a period of 0.5 to 
20 days. The high secondary masses result from mass ac- 
cretion during the second Roche-lobe overflow phase of 
the primary. 



3.2.3. Case C CE phase with a naked helium star or a 
white dwarf remnant 

Similar to formation channel 3, the initial orbital sepa- 
rations of the binaries following the sixth evolutionary 
channel are wide enough to avoid any type of Roche-lobe 
overflow until the primary reaches the AGB. However, the 
primary here fills its critical Roche lobe before the wind 
has a chance to reduce the primary's mass below 2 /3 that 
of its companion so that the resulting mass-transfer phase 
is now dynamically unstable. The binary emerges from 
the ensuing common-envelope phase as a WDMS binary 
consisting of a C/O or O/Ne/Mg white dwarf, the for- 
mer core of the AGB primary, and a low- to intermediate- 
mass main-sequence star. If the common-envelope phase 
takes place early on the AGB, the C/O or O/Ne/Mg core 
may retain a thin helium envelope which is subsequently 
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Fig. 8. As Fig. ^ but for formation channel 6. 



stripped away by a Wolf-Rayet type stellar wind. The 
main evolutionary phases occurring in the formation chan- 
nel are summarised schematically in Fig. [HI We will refer 
to this channel as formation channel 6. 

From the bottom panels of Figs. [21 and [21 it follows 
that the systems evolving through formation channel 6 
originate from binaries with initial orbital periods be- 
tween 60 and 6500 days, initial primary masses between 
1 Mq and 9 Mq, and initial secondary masses smaller than 
or equal to 8Mq. The lower limit on the orbital period 
range stems from the requirement that the binaries sur- 
vive the common-envelope phase, while the upper limit 
corresponds to the longest orbital period for which the 
primary fills its Roche lobe on the AGB. The limits on 
the primary and secondary masses arise for similar rea- 
sons as those in channel 3 (see Section ed. 3[) . 

The WDMS binaries forming through formation chan- 
nel 6 typically consist of a 0.5- 1.44 Mq C/O or O/Ne/Mg 
white dwarf and a main-sequence star with a mass up 
to 8Mq. The orbital period ranges from 0.05 to 2000 
days, where the long-period tail corresponds to systems 
in which the envelope mass of the primary is negligible 
in comparison to its core mass. A similar long-period tail 
was found by de Kool & Ritter (1993, see their Fig. 2a). 
Most WDMS binaries forming through this channel have 
Mwd « 0.5 - 0.6 M Q , M 2 £ 2 M Q , and P orb £ 20 days. 

3.3. Non-interacting systems 

The last formation channel is the most straightforward one 
as it represents the non-interacting WDMS binaries. The 
initial orbital periods of these systems are wide enough 
for the two stars to evolve in much the same way as they 
would if they were single. Their contribution to the popu- 
lation is therefore independent of the assumptions adopted 
for the treatment of mass transfer in semi-detached bi- 
naries, so that they may provide a convenient means to 
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renormalise our results for comparison with observations 
and with other authors. 

The two-dimensional distribution functions describing 
the population of WDMS binaries forming without in- 
teracting are displayed in the upper panels of Figs. El 
and [3| The initial orbital periods are typically longer than 
400 days, while the initial primary and secondary masses 
range from 1M to 9M and from 1 M© to 8M©, re- 
spectively. At the time of formation, most binaries have 
M WD w 0.5 - 0.6 M Q and M 2 & 2 M@. The final orbital 
periods may be substantially longer than the initial ones 
due to the action the stellar wind responsible for exposing 
the primary's core as a C/O or O/Ne/Mg white dwarf. 

3.4. The entire population 

In order to get an idea of the relative importance of the 
different formation channels, the formation space of all 
WDMS binaries forming through channels 1-7 is shown 
in Fig. [5] without renormalising the contributions of the 
different channels as in Figs. El and [31 The population is 
then clearly dominated by wide non-interacting systems 
with low-mass main-sequence stars evolving through for- 
mation channel 7. The second largest group consists of 
systems with low-mass main sequence stars and periods 
in the range from 1 to 10 days originating from channels 4 
and 6. The dominance of WDMS binaries with wide non- 
interacting progenitors furthermore implies that the ma- 
jority of the systems contain a C/O white dwarf with a 
mass around ~ 0.6M Q . Systems with white dwarf masses 
higher than ~ 0.8 M© or secondary masses higher than 
2 M© are relatively rare. This is further illustrated by the 
one-dimensional distribution functions for M\vd, M2, and 
Porb displayed in Fig. The bimodal nature of the or- 
bital period distribution is in good agreement with Fig. 5 
of de Kool & Ritter (1993). 

We will quantify the relative contributions of the dif- 
ferent formation channels in more detail in Section 

4. The effect of 7rlof and «ce 

The properties of the population of WDMS binaries 
forming through the evolutionary channels described in 
Section [31 may depend sensitively on the assumptions un- 
derlying the binary evolution calculations. The two main 
uncertainties applying to the discussed formation chan- 
nels are the amount of mass accreted by a normal-type 
star during episodes of stable Roche-lobe overflow and 
the treatment of the common-envelope phase resulting 
from dynamically unstable Roche-lobe overflow. In order 
to assess how these uncertainties affect the WDMS bi- 
nary population, we repeated the calculations presented 
in Section[3]for different mass-accretion parameters 7rlof 
and common-envelope ejection efficiencies acE- The val- 
ues of 7rlof and acE adopted in the different population 
synthesis models are listed in Tabled Our standard model 
(model A) corresponds to the assumptions described and 
used in Sections [3 and [31 



Table 1. Population synthesis model parameters. 
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The effects of 7rlof and «ce on the parameter space 
occupied by WDMS progenitors at the start of their evolu- 
tion and by WDMS binaries at the time of their formation 
are summarised in Table FJ1 As in Figs. [3 and [31 only bina- 
ries contributing to the present-day Galactic population 
are taken into account. In order to focus on the bulk of 
the systems forming through each formation channel, the 
mass and period ranges quoted are limited to bins con- 
taining at least ^1% of the channel's contribution to the 
WDMS binary population. In the case of model A, the ex- 
tent of the parameter space listed in the table is therefore 
somewhat more restricted than in Section [31 We also note 
that by simply quoting the minimum and maximum val- 
ues of Mi , Mi , and P or b we neglect the deviations of the 
high-density regions in the parameter space from a simple 
'rectangular' shape. 

Since the mass- accretion parameter 7rlof determines 
the amount of mass and thus the amount of orbital an- 
gular momentum lost from the system during episodes 
of dynamically stable Roche-lobe overflow, the parameter 
directly affects the evolution of WDMS binaries forming 
through formation channels 1-3. The parameter also in- 
directly affects the formation of WDMS binaries through 
channels 4-6 because it enters the stability criterion sep- 
arating dynamically stable from dynamically unstable 
Roche-lobe overflow systems (for details see, e.g., the ap- 
pendix in Kolb et al. 2001). The dependency of the crite- 
rion on 7rlof is such than an increase in 7rlof generally 
increases the critical mass ratio Q c = Ma onor / M &CCTCtoI 
separating dynamically stable (Q < Q c ) from dynamically 
unstable (Q > Q c ) Roche- lobe overflowing systems 3 . As 
we will see, the overall effect of the change in the stability 
criterion on the population is, however, rather small. 

From Table [3 it follows that the parameter space oc- 
cupied by the bulk of the WDMS binaries and their pro- 
genitors is fairly robust to changes in the mass-accretion 
parameter 7rlof- Channels 2 and 3 are the most sensitive 
to the value of the mass-accretion parameter. The most 
obvious overall effect of changing 7rlof is the decrease 
of the secondary mass M2 at the birth of the WDMS 

Note that the definition of the mass ratio as Q — 
Mdonor/Maccretor may differ from the definition of the initial 
mass ratio given by q = M2/M1. 
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Fig. 9. Distribution of present-day WDMS binaries in the (Mwd, P rb)- and (M2, P rb)-planes at the start of the 
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Fig. 10. White dwarf mass, secondary mass, and orbital period distributions of present-day WDMS binaries at the 
start of the WDMS binary phase. The distributions are normalised so that the integral over all systems found is equal 
to one. 



binaries with increasing degree of non-conservativeness. 
In addition, for a given initial primary mass M\, the in- 
crease of the critical mass ratio Q c with increasing values 
of 7rlof implies that in models G05 and Gl lower initial 
secondary masses become available for dynamically stable 
Roche-lobe overflow. Due to the shape of the initial mass 
function [Eq. (j^J] and the adopted initial mass ratio dis- 
tribution [Eq. J3J], this behaviour shifts the bulk of the 
systems evolving through formation channels 2 and 3 to- 
wards lower primary and secondary masses. An increase 
in 7rlof is usually also accompanied by a decrease of the 
minimum orbital period at the birth of the WDMS bi- 
naries because more angular momentum is lost from the 
system. The stellar masses and orbital periods resulting 
from fully conservative mass-transfer (model GO) are fur- 
thermore very close to those found in our standard model 
(model A). This concordance arises from the small values 



of 7rlof inferred from Eq. J2J) for binaries in which the 
thermal time scale of the accretor is not too much longer 
than the mass-transfer time scale of the donor. Finally, we 
note that the systems evolving through formation chan- 
nels 1-3 have typical orbital periods longer than 1 day, in 
agreement with the orbital period distribution of post- 
mode II WDMS binaries derived by de Kool & Ritter 
(1993). We recall that the terminology of mode II mass 
transfer was introduced by Webbink (1979) to indicate 
mass transfer from donor stars with radiative envelopes. 

The common-envelope ejection efficiency q;ce only af- 
fects the WDMS binaries forming through evolutionary 
channels 4-6. Since smaller values of cvqe require more 
orbital energy to expel the envelope from the system, the 
minimum orbital period required for a binary to survive 
the spiral-in process increases with decreasing values of 
acE- For a given donor star mass Mi, the donor star in 
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Table 2. The effect of the population synthesis model parameters on the masses and orbital periods of WDMS 
binary progenitors at the start of their evolution and of WDMS binaries at the time of their formation. For each 
formation channel, the mass and period ranges refer to bins containing at least ~1% of the channel's contribution to 
the WDMS binary population. The quoted ranges should therefore be considered as approximate ranges for the bulk 
of the systems and not as strict limits. For the construction of the table, only binaries contributing to the present-day 
Galactic population are taken into account. Details of the parameters adopted in the different population synthesis 
models are listed in Tabled 
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7RLOF = 1) 


1 


- 2 


< 2 


30 - 600 


0.3 


-0.5 


< 2 


0.1 


-20 


5 


CE1 


(a C E = 0.2) 


8 


- 11 


4-7 


30 - 100 


0.95 


- 1.4 


7-11 


2 


- 6 




A 


standard) 


5 


- 11 


2-9 


20 - 1000 


0.75 


- 1.44 


3-11 


0.5 


- 20 




CE8 


(oce = 5.0) 


3 


- 11 


< 7 


10 - 1000 


0.4- 


- 1.44 


< 9 


0.3 


-30 




GO 


7RLOF = 0) 


6 


- 11 


2-9 


20 - 1000 


0.8 - 


- 1.44 


3-11 


1 - 


- 20 




G05 


(7RLOF = 0.5) 


4 


- 11 


2-9 


20 - 1000 


0.7- 


- 1.44 


2-10 


0.5 


- 20 




Gl 


(7RLOF = 1) 


6 


- 11 


2-9 


200 - 1000 


0.8- 


- 1.44 


2-9 


1 - 


-20 


6 


CE1 


(oce = 0.2) 


1 


- 2 


< 3 


300 - 2000 


0.5 


-0.8 


< 3 


0.1 


- 10 




A 


standard) 


1 


- 5 


< 3 


100 - 2000 


0.5 


-0.9 


< 3 


0.1 


-60 




CE8 


(oce = 5.0) 


1 


- 5 


< 3 


30 - 2000 


0.5 


-0.9 


< 3 


0.2 


-300 




GO 


7RLOF = 0) 


1 


- 5 


< 3 


100 - 2000 


0.5 


-0.9 


< 3 


0.1 


-60 




G05 


(7rlof = 0.5) 


1 


- 5 


< 3 


100 - 2000 


0.5 


-0.9 


< 3 


0.1 


-60 




Gl 


7RLOF = 1) 


1 


- 5 


< 3 


100 - 2000 


0.5 


-0.9 


< 3 


0.1 


-30 


7 


A 


standard) 


1 


- 5 


< 2 


10 a - 10 8 


0.5 


-0.9 


< 2 


3 x 10 a 


- 5 x 10 s 



a binary surviving the common-envelope phase therefore 
has more time to evolve when cvce is small, so that the 
minimum mass of the nascent white dwarf also increases 
with decreasing values of cvce- Larger values of cxce, on 
the other hand, yield wider orbital separations at the end 
of the common-envelope phase, making it easier for short- 
period binaries to survive the spiral-in process. These ten- 
dencies are clearly seen in Table |2 when comparing the 
variations in the parameter space between models CE1 
(oce - 0.2), A (a CE = 1.0), and CE8 (a CE = 5.0). 

For a given a given orbital separation a; at the onset of 
the common-envelope phase and a given donor star mass 
Mi, the orbital separation at at the end of the common- 
envelope phase furthermore increases with increasing sec- 
ondary masses M2 ■ The mass of the secondary may there- 
fore be the deciding factor determining whether or not 
a binary close to the borderline separating merging from 
non-merging systems survives the common-envelope phase 



or not. From Eq. l[T|l. it follows that the minimum sec- 
ondary mass M2 required to survive the common-envelope 
phase increases with decreasing values of «ce- This be- 
haviour is responsible for the variations in the initial sec- 
ondary mass range of binaries evolving through formation 
channel 5. Since M1/M2 > Q c is a necessary requirement 
for the occurrence of a common-envelope phase initiated 
by the primary, the changes in the initial secondary mass 
range are accompanied by similar changes in the initial 
primary mass range. 

The influence of the donor star's mass Mi on the out- 
come of the common-envelope phase can be understood 
by noting that the radius of giant-type stars is predom- 
inantly determined by the mass of their core, and that 
the radius of the star at the onset of Roche-lobe overflow 
is approximately equal to the radius of their Roche lobe. 
For a given initial orbital separation a; and a given sec- 
ondary mass M2, the orbital separation af at the end of the 
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common-envelope phase therefore increases with decreas- 
ing values of M\. Consequently, the maximum primary 
mass Mi for which a system may survive the common- 
envelope phase and the maximum mass of the resulting 
white dwarf decrease with decreasing values of acE • This 
behaviour is responsible for the lower upper limit on the 
initial primary mass range of formation channel 6 when 
a CE = 0.2 (model CE1). 

Besides these general tendencies, the following addi- 
tional differences between the various models may be ob- 
served upon inspection of the individual formation chan- 
nels: 

• Channel 2. — In model GO, the minimum initial 
primary mass Mi increases from 2 Mq to 4 Mq . The in- 
crease is associated with the destabilising effect of smaller 
TRLOF-values on the mode of mass transfer from the pri- 
mary: since the main-sequence radius of a star increases 
with increasing mass, a lower mass star has more time 
to evolve and thus to develop a deep convective envelope 
before it fills its Roche-lobe than a higher mass star in 
a binary with the same orbital period. For 7rlof = 0, 
primaries with mass Mi <, 4M Q therefore lead to dy- 
namically unstable mass transfer which usually results in 
the merger of the primary's core with its main-sequence 
companion. We furthermore note that for model A a gap 
occurs in the initial primary mass range between 3 Mq 
and 4 Mq , in the white dwarf mass range at the time of 
formation between 0.4 M Q and 0.7 M , and in the sec- 
ondary mass range at the time of formation between 3 Mq 
and 6 Mq . This behaviour is caused by the dependency 
of 7rlof on the ratio of the mass-transfer time scale of 
the donor to the thermal time scale of the accretor [see 
Eq. As a consequence, some combinations of masses 
and orbital periods in model A are subjected to the same 
destabilising effect of small 7r.lof -values as described for 
model GO. Finally, we point that models A, G05 and Gl 
yield C/O white dwarfs with masses as low as ~ 0.3 Mq. 
Comparably low-mass C/O white dwarfs were also found 
by Iben & Tutukov (1985) and Han et al. (2000) on the 
basis of more detailed numerical calculations. 

• Channel 3. — In model Gl, the minimum initial pe- 
riod of WDMS binary progenitors is significantly shorter 
than in any of the other population synthesis models con- 
sidered. The reason for this is that in models A, GO, 
and G05, systems with initial orbital periods shorter than 
~ 600 days undergo a common-envelope phase instead of 
a dynamically stable Rochc-lobc overflow phase on the 
first giant branch, so that they evolve through formation 
channel 4 or 5 rather than through formation channel 3. 
When mass transfer becomes highly non-conservative, as 
in model Gl, the mass-transfer phase becomes dynami- 
cally stable and the primary can avoid losing its entire 
envelope until it reaches the AGB. 

• Channel 5. — In model G05, the minimum initial 
primary mass Mi is about 2 Mq lower than in mod- 
els GO and Gl. This behaviour is again associated with 
the stabilising effect of larger 7RLOF-values. Binaries with 
primary masses near 4 Mq do not survive the common- 



envelope phase in model GO, while they undergo a dy- 
namically stable mass-transfer phase in model Gl. In the 
latter case, the systems therefore evolve through formation 
channel 2 rather than through formation channel 5. When 
7rlof ~ 0.5 (model G05), mass transfer from primaries 
with Mi w 4 Mq initially takes place on the thermal time 
scale of the donor star, so that the dynamical instability 
is delayed until the donor reaches the red-giant branch. 
Contrary to the outcome of the common-envelope phase in 
model GO, the decrease of the donor star's mass during the 
initial thermal time scale mass-transfer phase here yields a 
primary mass at the onset of the common-envelope phase 
that allows the binary to survive the spiral-in process. For 
a somewhat related reason, the longest possible initial or- 
bital period decreases from ~ 1000 days to ~ 100 days in 
model CE1. As discussed in Section f3.2.2l systems with 
-Porb Ss 100 days start mass transfer on the thermal time 
scale of the donor star, while systems with P or b 100 
days go straight into the common-envelope phase. When 
O-CE = 0.2, the primary masses of the latter systems at 
the onset of the common-envelope phase are too high to 
avoid a merger during the spiral-in process, so that only 
the group with P rb Si 100 days survives. Finally, the in- 
crease of the lower limit on the initial orbital periods in 
model Gl is associated with the disappearance of bina- 
ries with Hertzsprung-gap donor stars (i.e. binaries with 
-Porb 55 200 days). For highly non-conservative mass trans- 
fer these binaries evolve through formation channel 2 in- 
stead of through formation channel 5. 

• Channel 6. — In model Gl, the upper limit on the 
orbital periods at the formation time of the WDMS bi- 
naries decreases with respect to models GO and G05. The 
stabilising effect of larger 7RLOF-values here implies that 
less mass needs to be lost via the stellar wind in order 
for mass transfer to be dynamically stable. Many of the 
longer period systems (i.e. the systems that have evolved 
furthest on the AGB and have thus lost the most mass 
in a stellar wind) therefore now evolve through formation 
channel 3 instead of formation channel 6. 

5. WDMS binary numbers and formation rates 

Similar to the variations in the parameter space discussed 
in the previous section, the assumptions adopted in the 
various population synthesis models listed in Tabled may 
affect the formation rates and the number of WDMS bi- 
naries currently populating the Galaxy. These variations 
do not apply to systems evolving through formation chan- 
nel 7 since their evolution is governed solely by mass-loss 
and mass-accretion from stellar winds. 

The total number of WDMS binaries currently popu- 
lating the Galaxy and the relative number of systems that 
formed through the interacting and the non-interacting 
formation channels are listed in Table for different ini- 
tial mass ratio or initial secondary mass distributions. The 
most striking effect is the strong decrease of the absolute 
number of systems for the initial mass ratio distribution 
n(q) oc q~°", for < q < 1. The decrease is related to 
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Table 3. Total number of WDMS binaries currently pop- 
ulating the Galaxy and relative contributions of the in- 
teracting and the non-interacting formation channels for 
different initial mass ratio or initial secondary mass dis- 
tributions. The absolute numbers may be converted into 
an approximate local space density of WDMS binaries by 
dividing them by 5 x 10 11 pc 3 (see Section l2~2*)l . 





Total Number 


Interacting Non-Interacting 


Model 


(Channel 1-7) 


(Channel 1-6) 


(Channel 7) 




n{q) 


= 1, < q < 1 




A 


1.6xlO y 


13.6% 


86.4% 


GO 


1.6xl0 9 


13.4% 


86.6% 


G025 


1.6xl0 9 


13.7% 


86.3& 


G05 


1.6xl0 9 


14.0% 


86.0% 


G075 


1.6xl0 9 


14.6% 


85.4% 


Gl 


1.7xl0 9 


15.5% 


84.5% 


CE1 


1.5xl0 9 


3.5% 


96.5% 


CE3 


1.6xl0 9 


9.8% 


90.2% 


CE6 


1.8xl0 9 


20.1% 


79.2% 


CE8 


1.8xl0 9 


23.4% 


76.6% 


n(q) <x q, < q < 1 


A 


1.3xl0 a 


15.6% 


84.4% 


GO 


1.3xl0 9 


15.4% 


84.6% 


G025 


1.3xl0 9 


15.8% 


84.2% 


G05 


1.3xl0 9 


16.2% 


83.8% 


G075 


1.3xl0 9 


17.2% 


82.8% 


Gl 


1.3xl0 9 


18.5% 


81.5% 


CE1 


1.2xl0 9 


4.9% 


95.1% 


CE3 


1.2xl0 9 


11.9% 


88.1% 


CE6 


1.4xl0 9 


22.4% 


77.6% 


CE8 


1.5xl0 9 


24.7% 


75.3% 




n{q) oc 


g-" ", < q < 1 




A 


6.6xlO Y 


10.7% 


89.3% 


GO 


6.6xl0 7 


10.6% 


89.4% 


G025 


6.6xl0 7 


10.7% 


89.3% 


G05 


6.7xl0 7 


10.9% 


89.1% 


G075 


6.7xl0 7 


11.2% 


88.8% 


Gl 


6.7xl0 7 


11.5% 


88.5% 


CE1 


6.0xl0 7 


2.0% 


98.0% 


CE3 


6.4xl0 7 


7.0% 


93.0% 


CE6 


7.3xl0 7 


18.6% 


81.4% 


CE8 


7.6xl0 7 


21.6% 


78.4% 


M 2 from IMF given by Eq. © 


A 


2.4xlO a 


10.3% 


89.7% 


GO 


2.4xl0 9 


10.1% 


89.9% 


G025 


2.4xl0 9 


10.3% 


89.7% 


G05 


2.4xl0 9 


10.4% 


89.6% 


G075 


2.4xl0 9 


10.6% 


89.4% 


Gl 


2.4xl0 9 


10.7% 


89.3% 


CE1 


2.2xl0 9 


1.8% 


98.2% 


CE3 


2.3xl0 9 


6.6% 


93.4% 


CE6 


2.6xl0 9 


18.2% 


81.8% 


CE8 


2.7xl0 9 


21.3% 


78.7% 



the larger number of systems undergoing dynamically un- 
stable mass transfer which ends in the merger of the two 
component stars. 

Although the absolute number of systems resulting 
from formation channel 7 is independent of the adopted 



population synthesis model parameters, the relative num- 
ber of systems changes due to variations in the number 
of systems forming through formation channels 1-6. The 
non-interacting systems generally account for 75% to 85% 
of the total population if the initial mass ratio is dis- 
tributed according to n(q) = 1 or n(q) oc q, for < q < 1. 
This relative number increases to 95% when cvce = 0.2 
(model CE1) due to the significant decrease of the number 
of systems surviving the common-envelope phase in forma- 
tion channels 4-6. In the cases where the initial mass ratio 
is distributed according to n(q) oc g _0 ", for < q < 1, 
or where the initial secondary mass is distributed in- 
dependently from the primary mass M\ according to the 
initial mass function given by Eq. @, the relative contri- 
bution of the non-interacting systems to the population of 
WDMS binaries may be even larger. 

In Table the contribution of WDMS binaries with 
interacting progenitors to the present-day Galactic popu- 
lation is further subdivided according to the followed for- 
mation channel. The absolute number of systems forming 
through formation channels 1-6 increases with increasing 
values of 7rlof, although the overall difference is rather 
small: the number increases by less than a factor of ~ 1.25 
between models GO (7rlof = 0) and Gl (7rlof = 1) for 
all initial mass ratio or initial secondary mass distributions 
considered. The effect of ojce on the absolute number of 
systems is larger: the number increases by about an order 
of magnitude between models CE1 (cvce = 0.2) and CE8 
— 5.0). A similar dependency on ctce was noted by 
Iben et al. (1997) for WDMS binaries with secondaries less 
massive than 0.3 M©. 

The relative contributions of formation channels 1- 
3 increase with increasing values of 7rlof due to 
the stabilising effect of non-conservative mass transfer. 
Correspondingly, the relative contributions of channels 4 
and 6 decrease with increasing values of 7rlof ■ The influ- 
ence of 7rlof is largest for formation channels 1-3 whose 
relative contributions increase by factors of ~ 5 — 15 be- 
tween models GO (7rlof = 0) and Gl (7rlof = 1). 
Increasing the common-envelope ejection efficiency acE, 
on the other hand, increases the relative contributions of 
channels 4 and 5. The increase is caused by the larger 
number of systems surviving the common-envelope phase 
resulting from case B Roche-lobe overflow events. The rel- 
ative contribution of formation channel 6, somewhat sur- 
prisingly, decreases with increasing values of «ce , but the 
absolute number of systems forming through this chan- 
nel still increases. The decrease of the relative number is 
therefore caused by the more prominent increase of the 
absolute number of systems forming through channels 4 
and 5. For similar reasons, the relative contributions of 
channels 1-3 also decrease with increasing values of «ce- 

The population of WDMS binaries with interacting 
progenitors is dominated by systems forming through for- 
mation channels 4 and 6. Depending on the adopted model 
parameters, they generally account for 80% to 95% of the 
population if the initial mass ratio is distributed according 
to n(q) = 1 or n{q) oc q, for < q < 1, or even for 90% 
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Table 4. Total number of present-day Galactic WDMS binaries that formed through formation channels 1-6 and 
relative contributions of each formation channel for different initial mass ratio or initial secondary mass distributions. 
The absolute numbers may be converted into an approximate local space density by dividing them by 5 x 10 11 pc 3 . 







Dynamically stable mass transfer 




Common- envelope phase 






Case B 


Case C 




Case B 


Case C 
i 






WD 


\, 

nHe star 


I 

nHe star 


/ 
WD 


\ 

nHe star 


I 

WD or nHe star 


Model 


Total Number 


(Channel 1) 


(Channel 2) 


(Channel 3) 


(Channel 4) (Channel 5) 


(Channel 6) 








n{q) 


= 1, < q < 1 








A 


2.2x10 s 


3.4% 


0.2% 


0.1% 


40.5% 


0.1% 


55.7% 


GO 


2.2xl0 8 


2.0% 


0.2% 


0.1% 


41.1% 


0.1% 


56.5% 


G025 


2.2xl0 8 


4.1% 


0.3% 


0.2% 


40.1% 


0.2% 


55.1% 


G05 


2.3xl0 8 


6.2% 


0.7% 


0.4% 


39.1% 


0.2% 


53.5% 


G075 


2.4xl0 8 


10.2% 


1.6% 


0.6% 


37.1% 


0.1% 


50.5% 


Gl 


2.6xl0 8 


14.8% 


2.8% 


1.5% 


34.4% 


0.1% 


46.5% 


CE1 


5.2xl0 7 


14.5% 


1.0% 


0.4% 


18.8% 


<0.1% 


65.2% 


CE3 


1.5xl0 8 


4.9% 


0.3% 


0.1% 


35.5% 


0.1% 


59.0% 


CE6 


3.7xl0 8 


2.0% 


0.1% 


0.1% 


47.1% 


0.4% 


50.2% 


CE8 


4.3x10 s 


1.8% 


0.1% 


0.1% 


49.7% 


0.8% 


47.6% 


n(q) oc q, < q < 1 


A 


2.0x10 s 


3.8% 


0.3% 


0.2% 


44.9% 


0.2% 


50.6% 


GO 


2.0x10 s 


2.7% 


0.3% 


0.2% 


45.4% 


0.2% 


51.2% 


G025 


2.1xl0 8 


5.2% 


0.5% 


0.4% 


44.2% 


0.2% 


49.6% 


G05 


2.1x10 s 


7.7% 


1.0% 


0.6% 


42.8% 


0.2% 


47.8% 


G075 


2.3xl0 8 


12.7% 


2.2% 


1.1% 


39.9% 


0.1% 


44.1% 


Gl 


2.5xl0 8 


18.6% 


3.6% 


2.6% 


36.0% 


0.1% 


39.3% 


CE1 


5.7xl0 7 


13.6% 


1.2% 


0.6% 


20.2% 


<0.1% 


64.3% 


CE3 


1.5xl0 8 


5.2% 


0.5% 


0.2% 


39.1% 


0.1% 


54.9% 


CE6 


3.2xl0 8 


2.4% 


0.2% 


0.1% 


52.6% 


0.5% 


44.1% 


CE8 


3.6x10 s 


2.2% 


0.2% 


0.1% 


55.6% 


0.9% 


41.1% 








n(q) oc 


g- u ", < q < 1 






A 


7.1xlO b 


2.3% 


0.1% 


<0.1% 


35.5% 


0.1% 


62.0% 


GO 


7.0xl0 6 


1.1% 


0.1% 


<0.1% 


35.9% 


0.1% 


62.7% 


G025 


7.1xl0 6 


2.5% 


0.2% 


0.1% 


35.4% 


0.1% 


61.8% 


G05 


7.2xl0 6 


3.8% 


0.4% 


0.1% 


34.9% 


0.1% 


60.8% 


G075 


7.5xl0 6 


6.0% 


0.8% 


0.2% 


33.9% 


<0.1% 


59.0% 


Gl 


7.7xl0 6 


8.6% 


1.6% 


0.6% 


32.6% 


<0.1% 


56.6% 


CE1 


1.2xl0 6 


13.3% 


0.7% 


0.2% 


17.9% 


<0.1% 


67.9% 


CE3 


4.5xl0 6 


3.6% 


0.2% 


0.1% 


32.0% 


<0.1% 


64.1% 


CE6 


1.4xl0 7 


1.2% 


0.1% 


<0.1% 


39.8% 


0.2% 


58.8% 


CE8 


1.6xl0 7 


1.0% 


0.1% 


<0.1% 


41.5% 


0.6% 


56.9% 


M 2 from IMF given by Eq. © 


A 


2.5xl0 8 


2.3% 


<0.1% 


<0.1% 


37.0% 


<0.1% 


60.7% 


GO 


2.4xl0 8 


1.0% 


<0.1% 


<0.1% 


37.5% 


<0.1% 


61.5% 


G025 


2.5x10 s 


2.3% 


<0.1% 


<0.1% 


37.0% 


<0.1% 


60.7% 


G05 


2.5x10 s 


3.5% 


0.1% 


0.1% 


36.5% 


<0.1% 


59.9% 


G075 


2.6x10 s 


5.4% 


0.2% 


0.1% 


35.7% 


<0.1% 


58.6% 


Gl 


2.6xl0 8 


6.9% 


0.4% 


0.3% 


35.0% 


<0.1% 


57.4% 


CE1 


4.1xl0 7 


14.0% 


0.1% 


<0.1% 


19.0% 


<0.1% 


66.9% 


CE3 


1.5xl0 8 


3.7% 


<0.1% 


<0.1% 


33.5% 


<0.1% 


62.8% 


CE6 


4.8x10 s 


1.2% 


<0.1% 


<0.1% 


40.9% 


<0.1% 


57.8% 


CE8 


5.8x10 s 


1.0% 


<0.1% 


<0.1% 


42.6% 


0.2% 


56.2% 



to 99% of the population if the initial mass ratio is dis- 
tributed according to n(q) oc q~ ", for < q < 1, or the 
initial secondary mass is distributed according to the ini- 
tial mass function given by Eq. @ . The short-period bina- 
ries resulting from channels 4 and 6 (see Table[2J together 



with the long-period systems forming from wide non- 
interacting progenitors are therefore responsible for the 
bimodal nature of the orbital period distribution shown 
in Fig. EH The smallest contributions to the population 
of WDMS binaries with interacting progenitors stem from 
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formation channels 2, 3 and 5. We note, however, that 
channels 2 and 5 are interesting candidates for the forma- 
tion of WD + early B-star binaries with long and short 
orbital periods, respectively (Fig.^J. We will address these 
channels in this context in a forthcoming paper (Willems 
et al., in preparation). 

In TableOO we list the present-day birthrates of WDMS 
binaries forming through each of the considered formation 
channels. The birthrate of systems forming through forma- 
tion channel 7 is independent of the adopted population 
synthesis model parameters and amounts to ~ 0.3 yr -1 , 
unless the initial mass ratio is distributed according to 
n(q) oc q~°", for < q < 1. In the latter case, the 
birthrate of wide non-interacting systems is of the order 
of 0.01 yr _1 . The birthrate derived for our standard model 
with n(q) = 1, for < q < 1, is in excellent agreement 
with the birthrate for C/O white dwarfs in wide binaries 
of 0.29 yr" 1 derived by Iben & Tutukov (1986b). 

The general dependencies of the birthrates of for- 
mation channels 1-6 on the adopted population synthe- 
sis model are similar to those of the absolute and rela- 
tive numbers of systems described above. The effect of 
changing 7rlof is largest for channels 2 and 3 where 
the birthrate increases by an order of magnitude between 
models GO (7rlof = 0) and Gl (7rlof = 1)- The impact 
of 7rlof on the other channels is typically smaller than a 
factor of ~ 3. Changes in cvce, on the other hand, predom- 
inantly affect formation channels 4 and 5. The birthrates 
of systems forming through these channels increase by 1 
and 3 orders of magnitude, respectively, between model 
CE1 (a CE = 0.2) and model CE8 (a C E = 5.0). In the case 
where the initial secondary mass is distributed according 
to the same initial mass function as the initial primary 
mass, the birthrate of channel 5 increases by 4 orders 
of magnitude. The birthrate of systems forming through 
formation channel 6 varies by less than a factor of ~ 5, 
while the birthrates of systems forming through formation 
channels 1-3 are not affected at all. For 7rlof 0.25 
and acE ^ 0.6, the total birthrate of all WDMS binaries 
forming through a common-envelope phase (channels 4- 
6) is furthermore about 10 times larger than the total 
birthrate of all WDMS binaries forming through a stable 
Roche- lobe overflow phase (channels 1-3). This is again in 
excellent agreement with de Kool & Ritter (1993). 

For conclusion, we consider the ratio of the number of 
systems containing a He white dwarf to the number of sys- 
tems containing a C/O or O/Ne/Mg white dwarf. This ra- 
tio is independent of the normalisation of the adopted star 
formation rate [Eq. ijfJJl] and of the fraction of stars in bi- 
naries. Due to the dominance of the wide non-interacting 
systems and the systems forming through the common- 
envelope channels 4 and 6, the ratio is fairly insensitive to 
changes in the mass-accretion parameter 7rlof- The role 
of the common-envelope ejection parameter acE, on the 
other hand, is much more prominent. This is illustrated 
in the left-hand panel of Fig. For cvce = 0.2 (model 
CE1) the ratio takes comparably small values for each of 
the initial mass ratio or initial secondary mass distribu- 



tions considered. Provided that a complete observational 
sample of WDMS binaries can be compiled and the na- 
ture of the white dwarf can be established, comparison of 
the observationally derived ratio of He white dwarf sys- 
tems to C/O or O/Nc/Mg white dwarf systems with the 
theoretically predicted values may therefore constrain the 
common-envelope ejection efficiency. For higher values of 
oce, the effect of changing cxce can not be distinguished 
from the effect of changing the adopted initial mass ratio 
or initial secondary mass distribution. As shown in the 
right-hand panel of Fig. 1111 similar dependency on 7rlof 
and acE is also found when the subset of short-period sys- 
tems resulting from the dominant channels 4 and 6 is con- 
sidered. This subset is potentially more interesting than 
the whole sample of WDMS binaries since shorter period 
systems are easier to detect observationally than longer 
period systems. 

6. Luminosity ratios and radial velocities 

The detection and identification of WDMS binaries gen- 
erally relies on the ability to detect the white dwarf's sig- 
nature in the composite binary spectrum and to deter- 
mine the components' orbital radial-velocity variations. 
Regardless of any additional selection effects, the de- 
tection of the white dwarf's spectral signature becomes 
increasingly difficult with decreasing ratio of the white 
dwarf's luminosity to the main-sequence star's luminosity, 
while the detection of radial-velocity variations becomes 
more difficult with increasing orbital periods. 

The luminosity ratio may be particularly stringent for 
binaries with older white dwarfs since the time during 
which a white dwarf is detectable is typically of the order 
of 10 s years or shorter (e.g. Iben et al. 1997). The distribu- 
tion of the luminosity ratios £wd/£mSj where Lwd is the 
luminosity of the white dwarf and Lms that of its main- 
sequence companion, for WDMS binaries forming through 
the different formation channels are displayed in Fig. ll2l in 
the case of population synthesis model A and the initial 
mass ratio distribution n(q) — 1. For the determination 
of the luminosity ratios, we adopted the expression for 
white dwarf cooling by Hurley et al. (2000) and neglected 
the relatively small increase of the secondary's luminosity 
during its evolution on the main sequence. The luminosity 
ratio distributions typically peak between 10~ 5 and 10~ 3 , 
except in the case of formation channel 4 where a much 
broader distribution is found which extends up to signif- 
icantly higher luminosity ratios. The figure furthermore 
shows that channels 4 and 6 contain a significant number 
of WDMS binaries with small luminosity ratios, so that it 
may actually be hard to compile a complete sample of sys- 
tems to constrain the common-envelope ejection efficiency 
acE- 

In Fig. 1131 we combine the luminosity ratio distribu- 
tions shown in Fig. 1121 with the relative contributions of 
the different formation channels to the WDMS binary 
population given in Tables E3 and 0] to determine the frac- 
tion of systems for wich the white dwarf to main-sequence 
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Table 5. Birthrates of Galactic WDMS binaries forming through formation channels 1-7 in the case of different initial 
mass ratio or initial secondary mass distributions. The rates may be converted into approximate local birthrates by 
dividing them by 5 x 10 11 pc 3 . 
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star luminosity ratio is larger than a given limiting ra- 
tio (£wD/'-£'Ms) cr ir The thick solid line represents the 
total for all systems forming through formation chan- 
nels 1-7, while the thin lines distinguish between the 
contributions of the individual formation channels. For 



(£wD/£Ms) cr it — 10~ 6 , the fractional contributions of 
the different channels correspond to those obtained from 
Tables 01 and 0] The wide non-interacting systems dom- 
inate up to (iwD/iMs) cr it ~ 10" 1 . For larger values of 




Fig. 11. Ratio of He white dwarf WDMS binaries to C/0 or O/Ne/Mg WDMS binaries for different initial mass 
ratio or initial secondary mass distributions. The left-hand panel shows the ratios as a function of the common- 
envelope ejection efficiency acE for WDMS binaries forming through all considered formation channels, while the 
right-hand panel shows the ratios for systems forming through formation channels 4 and 6 only. In the left-hand panel, 
the curves associated with n(q) oc q~°" and with Mi drawn independently from the same IMF as M\ are almost 
indistinguishable . 
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Fig. 12. Normalised distribution functions of white dwarf to main-sequence star luminosity ratios for formation chan- 
nels 1-7 in the case of population synthesis model A and the initial mass ratio distribution n(q) = 1. 



(Lwd / Lms) CI [ t i the systems forming through formation 
channel 4 become equally important. 

The amplitudes i^wD and Kms of the white dwarf 
and the main-sequence star's orbital radial-velocity vari- 
ations in the case of a typical orbital inclination i = 60° 
and under the assumption of a fixed Keplerian orbit are 
listed in Table [BJ This assumption breaks down for short- 
period systems in which magnetic braking and/or grav- 
itational radiation cause the components to spiral-in to- 



wards each other. However, we do not expect this to sig- 
nificantly affect our radial-velocity estimates, unless the 
orbital evolution time scale becomes shorter than or com- 
parable to the white dwarf cooling time. As present-day 
radial-velocity surveys easily reach accuracies down to 10 
kms -1 , all formation channels except channels 3 and 7 
produce WDMS binaries with readily detectable radial- 
velocity amplitudes. If the spectra of the component stars 
can be disentangled, the common-envelope scenarios fur- 
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Fig. 13. The fraction F of systems for wich the white dwarf to main-sequence star luminosity ratio is larger than 
the limiting ratio (£\yD/£Ms) cr j t j m ^ ne case of population synthesis model A and the initial mass ratio distribution 
n(q) = 1. The thick solid line represents the total for all systems forming through formation channels 1-7, while the 
thin lines distinguish between the contributions of the different formation channels. The inserts zoom in on the regions 
where the distribution functions are small. 



Table 6. Typical orbital radial-velocity amplitudes of WDMS binary components in the case of population synthesis 
model A and the initial mass ratio distribution n(q) — 1. The radial- velocity amplitudes correspond to a typical orbital 
inclination of 60°. 



Formation channel 



A v 



K 



MS 



1 (Stable case B RLOF -> WD) 

2 (Stable case B RLOF -> nHe star) 

3 (Stable case C RLOF -» nHe star) 

4 (Case BCE^ WD) 

5 (Case BCE^ nHe star) 

6 (Case CCE^ nHe star or WD) 

7 (Non-interacting systems) 



50 



100 kms" 
100 kms" 
5—15 kms" 
70 - 220 kms" 
170 - 220 kms" 
40 - 170 kms" 
1 - 10 kms" 



5-10 kms" 
5-10 kms" 
2 — 5 kms" 
30 - 140 km s" 
20 - 30 km s" 
20 - 140 kms" 
1 — 5 kms" 



thermore produce WDMS binaries in which both compo- 
nents may have observable radial-velocity variations. 



7. Concluding remarks 

We used the BiSEPS binary population synthesis code 
described by Willems & Kolb (2002) to study the popula- 
tion of detached white dwarf main-sequence star binaries 
forming through seven distinct evolutionary channels. In 
six of these, the white dwarf is formed through binary in- 
teractions resulting from Roche-lobe overflow. The chan- 
nels mainly differ in the stability and the outcome of the 
mass-transfer phase that eventually gives rise to the for- 
mation of the white dwarf. The six formation channels 
are depicted schematically in Figs. and The sev- 
enth evolutionary channel is characterised by the absence 
of any type of binary interactions other than mass and 
angular momentum exchange via a stellar wind. Binaries 
evolving through this channel are typically very wide, so 
that the white dwarf forms in much the same way as it 
would if the primary were a single star. 



Our results show that the wide non-interacting bina- 
ries generally comprise more than 75% of the total WDMS 
binary population. The remaining part of the population 
is dominated by systems undergoing a common-envelope 
phase when the white dwarf progenitor ascends the first gi- 
ant branch or the asymptotic giant branch. The total num- 
ber of systems currently populating the Galaxy, the rela- 
tive number of systems evolving through each formation 
channel, and the birthrates of systems forming through 
each channel are given in Tables [3HS1 f° r different differ- 
ent assumptions about the fate of the mass transferred 
during dynamically stable Roche-lobe overflow, different 
common-envelope ejection efficiencies, and different ini- 
tial mass ratio or secondary mass distributions. The total 
number of systems is of the order of ~ 10 9 when a flat 
initial mass ratio distribution n(q) = 1, for < q < 1, 
is considered. This number decreases to ~ 7 x 10 7 for an 
initial mass ratio distribution of the form n{q) ex q~ 99 , 
for < q < 1. An even stronger decrease is found for 
the birthrates of systems forming through the different 
formation channels. Besides the different assumptions re- 
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garding the treatment of dynamically stable and unstable 
mass transfer, the total number of WDMS binaries cur- 
rently populating the Galaxy also depends on the finite 
life time of the systems which, in the case of short-period 
systems, may be affected by the adopted magnetic braking 
prescription. A detailed analysis of the effects of different 
braking laws is, however, beyond the scope of this inves- 
tigation. 

The total number of Galactic WDMS binaries evolv- 
ing through a common-envelope phase is of the order 
of ~ 2 x 10 8 when a flat initial mass ratio distribution 
n(q) = 1 is considered, and of the order of ~ 7 x 10 6 when 
an initial mass ratio distribution of the form n(q) oc q~°" 
is considered. Under the assumption that the Galaxy has 
an effective volume of 5 x 10 11 pc 3 , these numbers corre- 
spond to approximate local space densities of 4x 10 -4 pc -3 
and 10 -5 pc -3 , respectively. The latter space density is 
comparable to the space density of 6 — 30 x 10~ 6 pc -3 de- 
rived by Schreiber & Gansicke (2003) from an observed 
sample of 30 post-common-envelope binaries with low- 
mass main-sequence secondaries. However, since the space 
densities derived by Schreiber & Gansicke (2003) are likely 
to be lower limits, the space density of 4 x 10~ 4 pc~ 3 is 
not necessarily in disagreement with observations. We also 
note that 12 of the 30 post-common-envelope binaries con- 
sidered by Schreiber & Gansicke (2003) contain an sdOB 
primary instead of a white dwarf. These sdOB + main- 
sequence star systems are not included in our simulated 
WDMS binary samples. 

We furthermore find that the ratio of the number 
of He white dwarf systems to the number of C/O or 
O/Ne/Mg white dwarf systems is weakly dependent on 
the amount of mass lost from the system during dynam- 
ically stable Roche-lobe overflow and strongly dependent 
on the common-envelope ejection efficiency. We there- 
fore propose that comparison of the observationally de- 
rived ratio from, e.g., large-scale planet-search mission as 
SuperWASP, COROT, and Kepler with the ratios pre- 
sented in Fig. 1111 may yield constraints on the common- 
envelope ejection efficiency, provided that the efficiency is 
low. For higher efficiencies, the role of the ejection effi- 
ciency parameter cannot be distinguished from the role 
of the initial mass ratio or initial secondary mass dis- 
tribution. We note, however, that the theoretically pre- 
dicted ratios are derived under the assumption of a con- 
stant envelope binding-energy parameter and a constant 
envelope-ejection parameter. In addition, the derivation 
of constraints on the common-envelope ejection efficiency 
requires the detection of WDMS binaries with luminosity 
ratios down to 10~ 5 , which poses a severe and not easily 
overcome observational challenge. 

This paper is to serve as a first comprehensive step 
in a study of more specific subclasses of WDMS binaries 
and binaries descending from them. In future investiga- 
tions we will address WDMS binaries in the context of 
pre-cataclysmic variables, double white dwarfs (Willems 
et al., in preparation), and white dwarf + B-star bina- 
ries (Willems et al., in preparation). Other potentially 



interesting applications are the study of type la super- 
novae progenitors through both the 'standard' evolution- 
ary channel and the recently proposed thermal time scale 
mass-transfer channel (King et al. 2003). 
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